Creation of spin 1/2 particles by an electric field in de Sitter space 
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Abstract 

In the present article we solve the Dirac equation in a de Sitter universe when 

a constant electric field is present. Using the Bogoliubov transformations, 

we compute the rate of spin 1/2 created particles by the electric field. We 

compare our results with the scalar case. We also analyze the behavior of 

the density of particles created in the limit H=0, when de Sitter background 

reduces to a flat space-time. 
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During the last two decades, quantum field theory in curved spaces has been extensively 
discussed, mainly in the study of particle creation in vicinity of black holes as well as in 
expanding cosmological universes. Here we are talking about a regime where the spacetime 
curvature remains significantly less than the Planck scale, and therefore one may usefully 
employ the semiclassical approximation in which the gravitational field is treated classically 
and we analyze quantized elementary particles propagating on this classical background. 
This is perhaps the simplest theory that one can use for discussing quantum effects in the 
presence of gravitational fields in the early universe. In this framework we are interested in 
discussing spin 1/2 particle creation by a constant electric field in a (1+1) de Sitter space. 
This problem has been preceded by a series of articles by Lotze |l],[| where the author 
analyzes the influence of interacting fields on the rate of particles produced by the gravity. 
In particular, he considers the electromagnetic interaction ie^f^A^ between the electron- 
positron field and the potential in expanding spatially flat Robertson- Walker Universes, 
and more recentlythe problem of pair production by an electric field in (l+l)-dimensional 
de Sitter space has been discussed for the case of scalar particles ||. Here, the idea was 
to establish a comparison with the results obtained after computing the rate of membranes 
created by an anti symmetric tensor field QJ5[] In the present article, as a natural extension 
of the results reported in Ref. || we solve the Dirac equation in (1+1) de Sitter universe 
when a constant electric field is present. After identifying the the "in" and "out" vacuum 
states, we compute, via the Bogoliubov transformations, the rate of particles created by the 
electric and gravitational fields. We compare our results with those obtained in the scalar 
case. 

The line element associated with a (1 + 1) de Sitter universe takes the simple form 



dg* = _ df 2 + e 2Ht dx 2 



(1) 



which, in terms of the conformal time r) 



— T?e reduces to 



1 



drj 2 + dx 2 ) 



(2) 
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Since we are interested in analyzing Dirac particle production by a constant electric field 
in a de Sitter background, in order to establish a comparison with the pure gravitational 
contribution, it is advisable to have as a previous result the rate of spin 1/2 pair production 
in de (1+1) Sitter universe. For this purpose we proceed to solve the Dirac equation in the 
(1+1) de Sitter space (0). The covariant generalization of the Dirac equation reads 

(l a (d a -r a )+m)V = (3) 

in the background field where 7° are the curved Dirac matrices satisfying the commutation 
relation {7°, 7^} = 2g a/3 , T a are the spinor connections, which for a diagonal metric as ([!]) 
reduces to 

Ta = -\g m TZ x sr, (4) 

where 

^ = ^(TY-7Y) (5) 

Choosing to work in the diagonal tetrad gauge, we have that the gamma matrices take the 
form 

7 ° = -#7/7°, 7 1 = -F777 1 (6) 

where the flat gamma matrices 7° satisfy {7°, 7^} = 2i] al3 . Then, from @, we obtain that 
the only non zero connection T a reads 

r, = -lyy (7) 

Since we are dealing with a (1+1) dimension space, we have that gamma matrices can be 
represented in terms of the Pauli matrices. In particular, we can make the following election: 

7° = i a \ f = a 1 (8) 

substituting the gamma matrices (0) as well as the spinor connection ([7]) into the Dirac 
equation @, and introducing the spinor $ related to \l/ as follows 



$ = V - 1/2 V (9) 

we obtain that the problem of solving Eq. @ reduces to find solutions of the coupled system 
of equations 

S + S*' + *-*»=° < io) 

where fc^ is the eigenvalue of the linear momentum operator —id x consequently, we can write 
down the spinor $ as $ = ^o(ri)e' lkxX Substituting ([11]) into ( |10D and vice-versa, we arrive 
at 



d 2 m 2 im 1 2 N 



(12) 



The solutions of the system of equations fll2[) can be expressed in terms of cylinder functions 
Z v (z) as follows 

$! = T] 1/2 Z v (k x T]), $ 2 = T^~l 2 Z v —\ (kx'n) (13) 

where the parameter i/ is given by the expression 

"=2 + lF (14) 

In order to analyze quantum effects in the background field (0) we have to define vacuum 
states in the asymptotic regions. The notion of positive and negative frequency states 
is associated with the existence of a timelike Killing vector. In expanding cosmological 
Universes that vector does not exist, and we have to appeal to another methods for obtaining 
a reasonable definition of the vacuum state. Among the different approaches we have that 
the so-called adiabatic method is perhaps the most popular. |||7| In our particular case, 
since the asymptotic behavior for scalar fields in the background field (0) was discussed 
in Ref. |§, we already know the asymptotic form that our in and out vacua should take. 



Nevertheless, a natural criterion for choosing the initial and final states is to look at the 
behavior of the quasiclassical modes described by the Hamilton- Jacobi equation 



g aP S, a S,p + m 2 = (15) 

which in the metric (0) can be written as a sum 

S = T( V ) + F(x) (16) 

and, after introducing the constant of separation k associated with the space variable x, we 
arrive at the first order differential equation: 

Then the quasiclassical behavior of the solution of the Dirac equation (|3|) in the metric (0) 
is 



, / exp(±i£; Jrf + ^r] ±j/m ' H 



(18) 

then, as rj — >• — oo (i — > oo) we find 

* -> ^o(^) exp(^ikrf) (19) 

and as r\ — > (£ — > oo) 

$ ^ ^(W^j-Wff^Wfl = $ W(^)-WS e T«/2 (20) 
if if 



where the upper and lower signs in ([191) and (pOj) are associated with positive and negative 



frequencies respectively. Then, looking at the asymptotic behavior of the Hankel functions, 

if < 2 ) fx) « J—e-**-™ Z 2 -/ 4 ) x — oo (21) 

V 7TX 

we have that the "in" state as rj — > — oo can be expressed in terms of H^(x) as 

*i = Tj 1/2 H<?\k ai v), (22) 



substituting (H) into (|32]) we readily find that $ 2 = -r] 1 / 2 ^ 2 }^^). The expression (\ZZ\ 



corresponds to the truly de Sitter invariant vacuum as it was pointed out by Bunch and 
Davies ||. From the definition of Hf?'(z) in terms of the Bessel function J v {z) 

Rf\z) = - J - v{z) (23) 

and the behavior for small values of z 

JJz) -> ^ -, (24) 

v ; 2T(z/+ 1)' V ; 

and as rj — > we have that in the vicinity of rj = the spinor component <J>i takes the form 
$i ^ ^ ((^T e *»T(l - i/) - 2T(i/ + (25) 



then, looking at the quasiclassical asymptotic behavior of the "out" states fl20|) and the 
linear decomposition between positive and negative frequency modes 

*&,k = a At,k + PkKut,k (26) 

we arrive at 

| ak |2 = e Wfl | p k |2 (27 ) 

following Mishima and Nakayama we can already establish that the spectrum of par- 
ticles created is a thermal one with temperature T = H/2tt From | | 2 + | | 2 = 1, one 
readily gets 

I P k ! 2= e 2,rm/tf + I ( 28 ) 

then, we obtain a thermal Fermi-Dirac distribution for | j3k | 2 , which for m » H gives as 
result that the number of particles per unit coordinate volume is ^ ~ e -2irm/H Notice 
that for the Dirac case, we do not need to impose any restriction on the relation between m 
and H in order to obtain well defined "out" states. 

Let us extend our problem to the case when a constant electric field minimally coupled 
to the spinor field is present. Then the Dirac equation can be written as 



(7 a («9 Q -r Q )+m-^e 7 a A Q )^ = 



(29) 



where the vector potential associated with a constant electric field E can be written in 
the comoving time r\ as follows 

A — 0, Aj — -§e»< = (30) 

substituting ( |30|) and T a given by (0) into ( |29"1) we find 

h+^-^K (3i) 



using the gamma's representation (H) we find that eq. (|3T|) yields 



o?77 -ffr? y 2 \ 



e£ 



$i = 



(33) 



the system of equations 



3T) can be written in the matrix form 







(34) 



and with the help of the similarity transformation T, acting on the Pauli matrices cr, as 
7V.T- 1 = a[ 

( 1 

-i 1 



T 



V 



;i-z<7i) 



(35) 



and the auxiliary spinor = rj 1 / 2 ®, we reduce the coupled system of equations (|32"D-(|3"3"D to 
the Whittaker equations 



d 2 1 ,ieE 11 +e 2 E 2 /H A + m 2 /H 2 + 1/4 X 
+ ± r - + s— ' — 



V^p 2 A H 2 T p 
having as a solution 







(36) 
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M A ,» = 9 1>2 = r, 1 ' 2 ^ = c^ff+^e-^Mifi - A + ~ 2p + 1, p) (37) 
where we have introduced the auxiliary variable p = 2ik x r], and the parameters p and A read 



o 



(38) 

where the upper and lower signs in fl38|) and fl3T|) correspond to $i and $2 respectively. A 
second set of solutions of the Eq. fl36| ) can be expressed in terms of the Whittaker functions 



W x ^{z) = z» +l l 2 e- z ' 2 U{ l - - A + p, 2p + l,z) 



(39) 



Using the quasiclassical asymptotae, we can choose positive and negative frequency modes. 
In fact, looking at the asymptotic behavior for large values of | z \ 



(40) 



and for small z 



(41) 



we obtain that ^ nk = cW\ tfl (z), where c is a constant of normalization, and $ 
M\ !fl (z). Then, using the relation || 

1 



M x Jz) = T(2p+l)e 



iir\ 



W.x^z) 



+ — — . 1 x exp 



r(p-A + |) r(p + A + | 



Z7r(A — jl — 



(42) 



3n it 

— < argz < -; 2p ^ -1, -2. 



from Eq. (f42|), and the relation 



(43) 



7T 



?/ sinh iry 



we arrive at 





2 


sinh 7r 


(Ia* 


-eE /H 2 ) 




-eE /H 2 ) 


1 A 


2 


sinh7r( p 


+eE /H 2 ) 





(44) 



(45) 



S 



the expression ( [45]) takes simpler form when | fi ± cEq/H 2 |>> 1 Taking into account the 



relation | a>k \ 2 + | (3k \ 
I Pk I 2 



1 we finally obtain the thermal distribution 



1 



3 (-f?[( e2E o+ m2 ^ 2 ) 1/2 - eE o]) 



(46) 



X _|_ e 2n\^\-2enE /H 2 

which shows that the temperature depends on the intensity of the electric field E . Also we 
have, as for the scalar case, that | (3k | 2 changes according the sign of k x . In fact, from fl32|) - 
(j33[), it becomes clear that changing the sign of k x is equivalent to replace e by — e in (|46|). 

It is interesting to analyze the form expression ([5]) takes in the limit H — > 0. In this 
case we have that the line element (|l|) reduces to a 2-d Minkowski (flat) metric. 
Since the exponent of | (3k | 2 in fl4*6"D can be rewritten as: 

-2tt eE 



-2tt r 



[e 2 E 2 n +m 2 H 2 ) l ' 2 -eEr 



and in the limit when H goes to zero, we have 

m 2 H 2 



lim [(1 + 



e 2 E 2 



a/2 



H 2 



1] 



1 | m2H \v> 

' e 2 E 2 > 



m 2 H 2 
2e 2 E 2 



we readily obtain that in the flat limit the density of particles created reduces to 



Pk |L= lim | (3, 



-■Km 2 1 (eEo) 



(47) 



(48) 



(49) 



which is the result obtained by Narozhny and Nikishov flll , [12|1 in analyzing the phenomenon 
of pair creation by a strong constant electric field when the transverse linear momentum p± 
is equated to zero. It is worth mentioning that expression (49) is proportional to the first 



term of the series obtained by Schwinger [pL3[] , giving the probability, per unit time, and per 
unit volume, that a pair is created by the constant electric field. 
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